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Abstract. We show that the study of the Calcium depletion is a valid and highly sensitive method for investigating 
the chemical and physical history of the very extended ionized nebulae seen around radio galaxies (EELR), massive 
' ellipticals and 'cooling flow' galaxies. By observing the near IR spectrum of nebular regions characterized by low 
— , excitation emission lines (LINER-like) , we can use the intensity of the [CaII]AA7291, 7324A doublet -relative to other 
£> ' lines, like Ha- to infer the amount of Calcium depletion onto dust grains. The presence of dust in these objects 
. -which does not necessarily result in a measurable level of extinction- would favour a 'galactic debris' rather than 
^""l ' a 'cooling flow' origin for the emitting gas. Before aplying such test to our data, we study four possible alternative 
. mechanisms to dust depletion and which could have explained the absence of the [Call] lines: a) ionization of Ca + from 
t—{ ' its metastable level, b) thermal ionization of Ca + , c) a high ionization parameter and/or a harder ionizing contiuum 
, than usually asummed and d) matter bounded models associated to a hard ionizing continuum. We show that none 
ON of these alternative mechanisms explain the absence of the [Call] lines, except possibly for the highly ionized EELR 
. where a high ionization parameter is required combined with a soft power law. We thus conclude that for the other 
low excitation emission regions (cooling flows, liners, low excitation EELR), the abscence of the Call lines must be 
q due to the depletion of Calcium onto dust grains. 

a: 

C/3 . 

> 

■ 1. Introduction 

The study of the interstellar medium (ISM) of external galaxies provides important information about the global 
kinematic (inflow, outflow) of such gas, its chemical composition and the implied star formation history, its mass 
distribution, etc. This gas forms a vital part of the record of the formation of the parent galaxy, and the evolutionary 
processes involved. 

How can this material be studied in details? One way is to have it illuminated or excited by a powerfull AGN, 
giving rise to the phenomenon of extended emission line region (hereafter EELR). The drawback of course is that it 
only allows us to look at a restricted class of galaxy. The large scale EELR phenomenon is observed in a majority of the 
most powerful radiogalaxies with EELR extending to radial distances of up to 100 kpc from the nuclei (Tadhunter 1986; 
Baum et al. 1988), much larger radii than the stellar population distribution of the parent galaxy. The morphologies 
and kinematics of such regions cover the full range from regular disc/ring systems to chaotic systems for which 
no pattern can be discerned. Their spectrum show strong emission lines, covering a wide range in ionization. It is 
generally accepted that the EELR are ionized by some mechanism connected with the nuclear activity, but there is no 
full concensus on the excitation mechanism. Some objects show evidences for an interaction between the radio jets, 
which transport energy to the outer radio lobes, and the gas in the outer region. Maybe this interaction is responsible 
of the excitation of the gas through some kind of shocks (Sutherland, Bicknell & Dopita 1993). For other objects, that 
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do not show evident spatial coincidence between the radio structures and the EELR, the excitation might be due to 
direct photoionization by the nuclear ionizing radiation field (Robinson ct al. 1987: hereafter RBFT87). A reduced 
scale version of the EELR is the one observed in many Seyfert galaxies (Haniff, Ward & Wilson 1988) where the ionized 
gas may extend up to a few kpc although it is brightest within the central 100-300pc. The EELR morphology tends 
to be conical (Wilson & Tsvetanov 1994). Its detailed observation is complicated by presence of the very luminous 
stellar background of the bulge. In normal bright ellipticals, quite weak extended nebulosities of low excitation is a 
common phenomenon (Buson ct al 1994, Goodfroij 1994). The very large scale gas around radio-galaxies and on which 
we focus here presents the advantage that the lines are observed against the sky background rather than against the 
bright parent galaxy bulge. 

The origin of the gas making up these large scale extended nebulosities remains unknown. Furthermore, the dis- 
tinction between these and the filamentary nebulae seen in some clusters and around some massive ellipticals and 
often identified with cooling flows, is unclear. We do know, however, that this gas in every case is chemically enriched 
as compared to primordial gas. Emission line analyses (RBFT87) show common element abundances to be within a 
factor of a few (<) of Solar and also to be rather uniform over all the objects observed. 

The two most likely explanations for the origin of the material are: 

a) debris from recent tidal interactions and mergers 

b) gas cooling from the hot (~ virial) phase which from X-ray observations (Forman, Jones & Tucker 1985) has 
been shown to exist around massive ellipticals and inside galaxy clusters 

The main arguments behind these explanations are the following: 

a) Heckman et al. (1986) showed that a large fraction of powerful radiogalaxies have morphological features - 
shells, tails, loops, etc- similar to those produced in numerical simulations of galaxy interactions (e.g., Toomre and 
Toomre 1972, Quinn 1984). This could indicate that the activity has been triggered either because fresh gas has been 
accreted from outside or because prexisting gas in the galaxy has been caused to collapse to the core as a result 
of the interaction. This interaction scenario is also supported by observations of a few nearby radio galaxies which, 
apart from morphological peculiarities, show large misalignments between the stellar and the gaseous rotation axes, 
indicative of an external origin for the gas. 

b) Hot X-ray coronae (T > 10 4 K) are a common feature of bright early-type galaxies. Within some critical radius, 
radiative cooling becomes important, leading to the development of the so-called 'cooling flow' hypothesis (Nulsen, 
Stewart & Fabian 1984, Thomas 1986, Thomas et al. 1986) Eventually, condensations or filaments could be formed, 
dense and cool enough to radiate detectable optical emission lines. Most powerful radio galaxies are too distant for the 
characteristic X-ray emission to be currently detectable, but it is quite plausible that the EELR gas has condensed out 
of a surrounding cooling flow. Any discovery of their EELR as consisting of extensive optical filaments would argue in 
favour of the cooling flow hypothesis. 

Distinguishing between these alternative hypotheses has been attempted using gas kinematic measurements (Tad- 
hunter, Fosbury & Quinn 1989) which show that the radio galaxy EELR generally have a high specific angular 
momentum which is difficult to reconcile with the cooling flow picture. An alternative approach is to look for the 
presence of dust associated with the EELR gas. If the gas has cooled directly from a hot phase, there will have been 
no opportunity for dust to form, according to the standard cooling flow theory. Any dust introduced from galaxies 
into the hot intraclustcr medium will be rapidly destroyed (Draine & Salpetcr 1979). We would need a mechanism 
to produce this dust once the gas has cooled down (Fabian, Johnstone and Daines 1994). If, on the other hand, the 
material has fallen in during a merger, the dust/gas ratio is expected to have a value appropriate to the gas chemical 
composition found in normal galaxies. 

Determining the presence or absence of dust is important because: 

-of its relevance in deriving chemical composition which takes into account the effects of depletion of metals unto 
dust grains 

-of the implication for the star formation history and the ISM evolution: when does dust form? (and so stars?) 

-of the effects of dust on the apparent morphology of continuum and line features: pure absorption (reddening) 
and/or scattering (blueing/polarization) 

We discuss below how the absence of the forbidden [Call] lines can be used to infer the presence of dust mixed 
in with the emission gas. We assess in detail all the most plausible alternative explanations to that of internal dust 
for explaining the absence of [Call] lines. As no acceptable alternative solution is found, we conclude in favour of 
the validity of the method initially proposed by Ferland (1992). We however adapt and optimize the [Call] dust 
detection technique to the context of the EELR studies in which we are involved. The observational results and their 
interpretation will be presented in a subsequent paper. 
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2. Outline of the method 

To investigate how the forbidden Calcium doublet of [Call] in the infrared is affected by the physical condition 
encountered in photoionized plasma, we outline first the computer code which we used and then proceed to illustrate 
how the [Call] lines might be used to infer dust in the ISM of galaxies and thus our interest in securing these conclusions 
by closing up the possibility of alternative interpretations. 

2.1. The photoionization code MAPPINGS. 

To compute the emission lines used in our study, we have used the multipurpose photoionization-shock code MAP- 
PINGS (c.f., Binette, Dopita & Tuothy 1985). The Ca+ ion is considered as a five level atom. Its structure is shown 
in Fig. 1. 

One interesting aspect of the updated (Binette et al 1993a, b) code MAPPINGS is that the effect of dust scattering 
on the line transfer is explicitly solved using the numerical solution of Bruzual, Magris & Calvet (1988) as described in 
Appendix A of Binette et al 1993b. Other effects of dust on the ionization structure as well as on the thermal balance 
of the plasma are considered followoing the prescriptions of Baldwin et al. (1991 see Appendix C). The dust grain 
charge is calculated self-consistently and the formula describing the photoelectron energy distribution and the yield 
are from Drainc (1978) but with a cap of 0.2 for the yield at high photon energies. 

We now proceed to define the parameters employed in the calculations. Most of these are derived from our obser- 
vational knowledge of EELR although we recognize that it is incomplete and can be biased by particular diagnostic 
tools which are employed. 

2.1.1. Metallicity 

For definitcness, we adopt a set of solar abundances (Anders & Grevesse 1989) for the trace elements, in line with 
the finding of RBFT87 who indicated values for the EELR not radically different from solar. The solar abundance of 
{not depleted) atomic Calcium is 2.2 x 10~ 6 by number relative to hydrogen. Depletion into interstellar dust grains 
is known to reduce this abundance in the local ISM by a factor of <~ 5000 (Whittet 1992) so that if depletion was 
taking place in the partially ionized EELR plasma, even a small dust-to-gas ratio might be sufficient to eliminate any 
detectable trace of atomic Calcium. 

2.1.2. Geometry 

It's been accepted for a long time that the NLR and the EELR are formed by individual clouds that are ionized 
by the central source. In most models presented here, each emitting cloud is considered as a radiation bounded slab 
(optically thick to the Lyman continuum) which comprises: (1) a fully ionized region, closer to the illuminated face 
and responsible for the high ionization lines; (2) a partially ionized zone (PIZ), where the low ionization lines like 
Call arc emitted. The boundary of the photoexcited regions (1 + 2) is defined as the depth where the following two 
conditions arc simultaneously satisfied: a) The unabsorbed ionizing flux is < 1 % of the impinging flux and b) the 
ionized fraction n H + /tih < 1 %. Matter bounded clouds will be considered in Sect. 3.4. As we want to concentrate on 
the depiction phenomena, wc will simplify the calculations by not considering any neutral region (3, see Fig. 2) beyond 
the PIZ which would contain dust, the effect of which would be to cause additional extinction for an hypothetical and 
unfavorably placed observer looking from the back of the slab. We have verified that the results reached here using 
line ratios of similar wavelengths are not altered by the presence or not of this neutral absorbing zone. 

2.1.3. Gas density 

We adopt a representative density for the EELR clouds of 300 cm~ 3 . Typical electron densities values measured vary 
between a few tens (or less) to a few hundreds. The densities derived from forbidden line ratios might not apply 
to every subregion, but they are sufficietly low to consider the low density limit as a generally valid and very good 
approximation to the physical conditions affecting the Call lines of large scale EELR. In the low density regime, the 
effects of density variations on the emission line ratios which are considered in this paper are very small by comparison 
to the effects of other parameters like the ionization parameter or the hardness of the ionizing continuum. 

The calculations consider the gas pressure to be constant within the cloud (isobaric models) with the density 
behaviour modulated with depth into the cloud by the behaviour of the temperature and by the ionization fraction of 
the gas. 
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2.1.4. The ionizing continuum and the ionization parameter U 

We implicitely assume that the dominant ionization mechanism of the ionized gas in radio galaxies is photoionization. 
It has been shown by RBFT87 that a hard continuum extending well down into the soft X-ray region, despite some 
discrepancies with the observed spectra, can be considered to reproduce generally well the measured line ratios. This 
continuum can take its source in the active nucleus or be locally generated (e.g., fast shocks: Bincttc, Dopita, Tuohy 
1985). We have considered power law (PL) distributions of various values of a (F„ cx v +a ) as well as hot blackbodies 
(BB) of temperature T bh rss 10 5 K in order to study the effects that hardness has on our conclusions. 

The ionization parameter, a measure of the excitation level of the ionized gas, is defined as the quotient between 
the density of impinging ionizing photons and the density of the gas cloud: 

x t Ionizing energy flux J L v dv jhv 

cntf 4:ivr 2 cnH 

where L v is the monochromatic luminosity of the source, r the distance of the cloud to the continuum source, nn the 
density of the gas and c the speed of the light. 

We find that the parameters having the strongest effect on the line spectrum are the ionization parameter U and 
the mean ionizing photon energy {i.e., hardness of the continuum). 

2.2. How to detect dust 

How can we detect dust within the gas associated with EELR? We here summarize different techniques and compare 
their sensitivities. 

a) Reddcnning: given the nature of the gas distribution and the fact that it is ionized externally (unlike HH 
regions which are internally excited), the extinction may not be necessarily sufficient high to be easily detected using 
optical observations. Indeed, the EELR ratios studied by RBFT87 show little or no reddening. Furthermore a small 
enhancement of the Balmer decrement over recombination case B might be interpreted as resulting from collisional 
excitation rather than from reddening. 

b) Scattering: polarization measurements of high redshift radio galaxies (di Serego Alighieri et al. 1989, 1992, 
Januzi & Elston 1991, Tadhunter et al. 1992) show conclusive evidence for scattered nuclear light over large volumes 
and, altough there are reasons to believe that the scattering medium is dust, it is difficult to rule out entirely Thomson 
scattering by hot electrons. Detailed studies of the low redshift galaxy PKS2152-69 do, however, show polarized 
continuum radiation from highly excited extranuclear gas cloud with an energy distribution which is so blue that it 
must arise from dust scattering (Fosbury et al. 1990). 

c) Infra-red thermal dust emission: the IRAS satellite has shown that many galaxies radiate significant fractions 
of their energy in the far infrared sprectral region. Significant masses of dust at temperatures of around 40K are 
responsible of this radiation at wavelengths of 60^m and beyond. In many cases the FIR spectral energy distributions 
is still rising at 100/im, out of the spectral range detectable by IRAS. The cool dust can only be detected at milimctrc 
and submilimetre wavelengths. A strong limitation of this technique is the very poor spatial resolution of the IRAS 
satellite. Groundbased studies of the far IR emission of galaxies in the sub- mm range also exist(e.g. Clements, Andreani 
& Chase 1993), but still with poor spatial resolution. Although it is possible to infer masses and temperatures of the 
warm and cool dust components (dependent on models), the spatial distribution of the dust is not known. 

d) Indirect effects on the line spectrum. There are several ways this can happen: effects of dust on the gas tempera- 
ture (photoionization of dust grains may raise the temperature of the plasma) , effects of dust on the ionisation structure 
(dust grains selectively absorb ionising photons of lower energies), and influence on apparent chemical composition 
via the depletion of refractory elements onto dust grains. The first effect does not provide a unique interpretation for 
the unusually high temperatures seen in some EELR (Tadhunter, Robinson & Morganti 1989) while the second effect 
cannot be discriminated against reliably since even dust-free photoionization models are still too uncertain to be used 
as absolute reference point. For these reasons, the last effect is the only clearly promising one and is looked into details 
below. 

Calcium is very sensitive to the presence of dust as it is always found to be depleted in the interstellar medium (e.g., 
Crinklaw, Federman & Joseph 1994). Photoionization calculations appropriate to LINERS -a hard ionizing spectrum 
with a relatively low ionization parameter- invariably predict the [CaII]AA7291, 7324A doublet (4s 2 S-3ePD) to be very 
strong (Ferland 1993). These two forbidden lines (hereafter Fl= A729lA F2= A7324A ) have a high critical density 
~ 10 6 cm~ 3 . The latter line, 7324A is the weakest of the doublet and is furthermore blended with the [OIIJA7325A 
multiplet. The other line, A729lA lies some 30A shortward of [Oil] and is therefore straighforward to isolate given 
reasonable spectral resolution. The fact that any of these doublet lines are generally not seen in LINERS but are so in 
some novae when their envelopes have reached the appropiate ionization level can be interpreted as evidence of Calcium 
depletion onto dust grains in the former objects. Since EELR which are the subject of our investigation are often seen 
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in the ionization parameter region of the line ratio diagnostic diagrams occupied by LINERS (10~ 4 <U< 10~ 3 ), we 
can similarly use the [Call] doublet measurements to infer whether or not there is depletion taking place in EELRs 
and thus conclude whether dust is also mixed with the gas as is thought to be the case in LINERs. 

The above arguments have already been used for several objects with 'cooling flow' filaments by Donahue & Voit 
(1993) to infer the presence of dust mixed with the ionized gas. Ferland (1993) has shown the great sensitivity of this 
method to the presence of dust under NLR conditions. We now show it to be also the case under EELR conditions. We 
present two diagnostic diagrams in Fig. 3, the ratio [Call] /Ha (7291/6563) and the ratio [CaII]/[ArIII] (7291/7135) 
as a function of the ionization parameter U . The variable parameter distinguishing the three different sequences in U 
is /z, the dust-to-gas ratio of the plasma expressed in units of the solar neighborhood dust-to-gas ratio. The dramatic 
difference in line ratios between the grain depleted Ca/H (/z > 0) and the undepleted case (/z = 0), shows the sensitivity 
of this method, particularly for low values of U . 

3. Possible alternative explanations to depletion. 

Before we carry the conclusions of the current analysis to the interpretation of our observations (Villar-Martm & Binette 
1995), we report first on our effort in investigating other possible alternative mechanisms to dust depiction. If the warm 
Ca + region predicted by standard models does exist, then Calcium depletion becomes the only reasonable explanation 
for the non detection of the doublet lines. What we consider in this section is the possible NON EXISTENCE of the 
emitting [Call] region by investigating different mechanisms which could eliminate it. During our investigation, we 
require however that successfull models do not result in important discrepancies with other observed line ratios. The 
mechanisms we have considered to eliminate the [Call] region are 

- Ionization of Ca + * by Lya and soft continuum photons from the metastable level of Ca + 

- Thermal (collisional) ionization of Ca + 

- Photoionization with a much harder continuum or a much higher U than usually asummed 

3.1. Ionization by Lya and soft continuum photons. 

Wyse (1941) proposed that the ionization of Call from the metastable level by Lya photons, could explain the fact 
that the IR lines of Call at 8498, 8542 and 8662A appear in emission near the maximum phase of Me variables, whereas 
the H and K lines only occur in absorption. Trapped Lya photons could also play a part in ionizing metastable Ca + * 
as suggested by Wallerstein et al. (1986). 

We investigate here if this process is important under the conditions found in EELR clouds. In order to do this, we 
add two terms to the ionization equilibrium equation of Call. One which considers photoionization of excited Ca + * 
by the impinging UV continuum. The other is photoionization of excited Ca + * by the nebular Lya photons. The 
statistical equilibrium equations give the relative population of the mestastable level, which turns out to be, under 
EELR conditions, 2ai ^ 10~ 7 , being n 4s the density of Ca + ions in the ground level. With such a negligible population, 
the density of ionizing photons must be very high to increase the ionization rate to a non negligible level as compared 
to the ground state ionization rate. A simple estimate presented in Appendix demonstrates this level to be out of 
reach. 

In summary, the very diluted radiation fields and the low densities appropriate to the EELR implies an extremely 
small population for the excited levels which prevents the ionization of Ca + * by Lya and soft continuum photons from 
being of any significance. 

3.2. Thermal ionization of Ca + . 

We investigate here the possibility of collisional ionization by thermal electrons of Ca + to Ca ++ , a process which 
is important when the electronic temperature becomes higher than 20000K (Jordan 1969). In order to establish a 
comparison in U, we have considered two extreme cases in our calculations, \ogU = —4 & — 2. To illustrate how a 
much harder continuum will result in much higher gas temperatures, we also use two different PL of index a — 1.4 
and —0.4. Note that such a hard continuum as a = —0.4 is probably quite unrealistic. However, our intention here is 
simply to test whether very high temperatures can be achieved with photoionization models and specifically near the 
Ca+ region. The results are shown in Fig. 4 as a function of depth in the photoionized slab. Of the eight plots, the 
four upper ones correspond to a = —1.4 while the four at the bottom to a = —0.4. The four plots on the right, have 
logU= —2, and the four on the left, logU= —4. Two plots therefore are shown for each pair of [U, a] values: one is the 
the temperature T4 in units of 10000K and the one immediately underneath is the intensity of Fl (erg.s .cm < ), 
both as a function of the depth in units of 10 20 cto into the slab. These plots allow us to see the correspondig electron 
temperature to the position where the bulk of the [CaII]Fl emission takes place. 
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For the traditional PL of index a = —1.4, the electron temperatures are not anywhere near high enough for the 
process of thermal ionization to be relevant. Harder continua and increasing U values do produce higher temperatures, 
but, unless U are unrealistically high (even logU= —2 is not enough), the gas is never hot enough. Thus, realistic 
photoionization models are not able to heat the gas sufficiently to thermally ionize Ca + . We might conjecture that 
there could be an additional heating source like shocks which could raise the temperature of the gas. This would be 
an interesting point for further investigation. 

3.3. Effects on Fl of varying U and the continuum hardness. 

We now investigate how the Fl line might become undetectable by simply varying arbitrarily the ionization para- 
meter or the continuum hardness. Fig. 5 shows six diagnostic diagrams with the absciss always representing the line 
ratio [OI]A6300/[OIII]A5007. The [OI]/[OIII] ratio monotonically increases with decreasing gas excitation (i.e., with 
decreasing U) and is therefore a good measure of the excitation level of the gas. Each diagram shows in ordinate a 
different line ratio which can be related to a given gas property. [01] A6300/Ha, for instance, might measure the hardness 
of the continuum. In the last diagram, the ordinate corresponds to the quotient F1/[OI]A6300. The three sequences 
of models shown in each diagram differ by the slope of the power law which takes on the values of a = —1.4,-1 and 
—0.4. The values of Log U covered by each curve is in the range —4 to —1. Our aim is to look for models which can 
decrease the Fl intensity below the detection limit. Let's look at how we might define a practical detection limit. The 
open squares in the diagrams of Fig. 5 represent line ratios measured by RBFT87 in several EELR. The faintest line 
they measure is typically HcIIA4686. The mean ratio of HeIIA4686/[OI]A6300 observed is 10~ 0,4 = 0.4 for the large 
scale nebulosities. We establish our 'artificial' detection limit in the following way: since Hell is one of the weakest 
line successfully measured by RBFT87, we will assume that any line fainter than 0.4 below the [OI]A6300 flux is not 
detectable. In the last diagram, the region where Fl falls below this detection limit is shown by a dash line grid. Any 
model found in this area is deemed successful in explaining the non-detection of Fl without requiring depletion. 

We see that for the standard PL (a = —1.4), only models with high U (log U> —2) decrease log(Fl/[OI]) below 
-0.4. These models, as we can see in the diagrams, would therefore be valid only for the high excitation EELR, but not 
for LINERs, cooling flow filaments, or EELR of low and intermediate excitation. On the other hand, increasing the 
hardness of the continuum (flatter power laws), helps Fl to get fainter with respect to [OIJA6300, but the discrepancies 
with observed line ratios in other diagrams become totally unacceptable (see top two diagrams of Fig. 5). 

It is interesting to compare a BB sequence (1.2 x 10 5 K) with the canonical PL sequence a = —1.4. We see that 
both ionizing continua reproduce rather well the observed line ratios as was earlier shown by RBFT87. From these 
line ratios alone, there are no reasons to favour power laws over hot blackbodies. A similar conclusion was reached by 
Bincttc, Robinson and Courvoisier (1988) for the mean NLR spectrum of Seyferts. 

A BB produces a much stronger Fl compared to [OI]A6300 than any of the power laws considered. One reason for 
this is that the fraction 0°/0 + in the PIZ is completely controlled by the charge exchange reactions of 0° and + 
with H + and H°, respectively, and not by direct photoionization of 0°. This is not so for Ca+/Ca° which is free to 
respond to the different amount of hard photons (the only one to make it to the PIZ) available in a PL or a BB. We 
conclude from the last diagnostic diagram that the BB models could not explain the abscence of Fl from the observed 
spectra without invoking depiction. 

3.4- Effects on Fl of truncating clouds 

We showed in the previous section that a harder continuum can potentially bring Fl under the detection limit but 
result in important discrepancies with the low excitation lines. The reason is that the partially ionized zone (PIZ) 
where most of the low excitation lines are generated gets larger and larger with increasing hardness of the continuum. 
If the clouds were truncated, the smaller PIZ would generate weaker low ionization lines, thus improving the overall 
fit. 

We have investigated models with a — —0.4 which were truncated at a depth which satisfies a given criterion based 
on a specific line ratio. This has been done in two ways: 

1) The criterion in this case is to truncate the calculations when OIII/H/3 has reached a value of 10 which is the typical 
ratio for the high excitation EELR. 

The sequence of models shown in Fig. 6 are separated by a factor of 0.14dex in U. Altough there are still discrep- 
ancies, there is a notable improvement compared to the radiation bounded models of Fig. 5. The predicted line ratios 
are now located closer to the observed data (same scale as in Fig. 5). We conclude that models with a hard ionizing 
continuum must be matter bounded in order to fit acceptably most observed ratios. 
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What happens now with Fl/[OI]? In the last diagram of Fig 6, we see that all these models produce Fl above the 
detection limits and cannot therefore explain the non detection of Calcium. 

2) The second method consists in imposing that the models produce Fl under the chosen detection limit (F1/[OI]=0.4) 
and check if such models agree with the position of the observed line ratios in the rest of the diagrams. 

The models in the sequence which satisfy this criterion are found in the range 3.7 x 10~ 3 < U < 2.7 x 10~ 2 . They are 
represented in Fig. 6 as open triangles connected by a solid line. As we see in the top left diagram, the ratio [OIIIJ/H/3 
is not any more defining a simple trend with excitation (which is represented by the ratio [OI]/[OIII]). Furthermore 
[01] /Ha remains a discrepant ratio as in Fig. 5. So although it is in principle possible to satisfy the F1/[OI]=0.4 
criterion with a hard power law, the truncation must be done at a specific yet ad hoc depth and furthermore the 
previous trend of excitation with U has disappeared. 

Without rejecting the possibility of a more complicate mixture of matter and radiation bounded models, we believe 
that simply truncating clouds does not convincingly solve the problem of the weakness of the Call doublet and dust 
depletion remains the most likely interpretation. 

It is interesting to note that truncated clouds adjust better the HeII/H/3 ratio (bottom left diagram) as proposed 
before by Morganti et al. (1991) and Viegas and Prieto (1992). 

4. Conclusions 

This work is based on the method proposed by Ferland (1993) to investigate the presence of dust mixed with the gas of 
the Narrow Line Region of active galaxies. Because photoionization models predict remarkably strong forbidden lines 
[CaII]AA7291,7324A assuming reasonable abundances of atomic Ca, the basic idea is to infer a systematic depletion 
of Calcium onto dust grains whenever the infrared [Call] lines are observed very weak or undetected. This test of the 
dust content was applied to cooling flow filaments by Donahue & Voit (1993) who concluded on the presence of dust. 

We have shown here that this sensitive method is also applicable to the conditions found in the EELR of radiogalax- 
ies. In order to make more secure any inference about the presence of dust based on [Call] lines, we have investigated 
alternative explanations for their absence: ionization of Ca + * to Ca ++ by Lya photons and soft continuum photons 
from the metastable level of Ca + , thermal ionization of Ca + , ionization of Ca + due to either a very high U value 
(ionization bounded case) or to a hard continuum (with truncated clouds). Except for the highly excited EELR which 
might not possess any Ca + region due to their high ionization level, the results are negative: none of the alternative 
mechanisms or models studied can explain the absence of the [Call] lines without dust depletion. 

Our conclusion is that the dust content test appears generally valid for the EELR of radiogalaxies (unless the 
excitation level of the gas is extremely high). This will allow us to make important conclusions about the origin of 
such gas, discriminating between galactic debris and the standard cooling flow theory 

A BB ionizing continuum characterized by a temperature of 1.2 x 10 5 K can reproduce the observed line ratios of 
the EELR at least as well as a PL of index a = — 1.4. On the other hand, the F1/[0I] from a BB is higher than that 
of a PL so the case in favour of depletion is even stronger. 

In a follow up paper, we will present long slit spectra of EELR, cooling flow filaments and Seyfert 2 NLR, all taken 
in the region of the [Call] doublet. The goal will be to apply the test of the Calcium depiction described above in 
order to conclude whether or not the gas in these nebulosities is mixed with dust. This will be our starting point for 
deciphering the origin of the emitting gas. 
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A. Appendix 

By comparing the estimated photoionization rates from the excited level 3d of Ca + *, H^d, to that from the ground 
state 4s of Ca + , H^, we show that photoionization of excited Ca + * is a negligible process. 

A fundamental parameter which determines the importance of the ionization rate from the metastable level is its 
relative population with respect the ground state: ^ai, being m s the density of Ca + ions in the ground level and 
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nsd the density of Ca + ions in the metastable level, 3d. To evaluate this fraction, we have solved analytically the 
statistical equilibrium equations. To simplify the calculations we have considered a three level atom, reducing the 
two 4p sublevels (see Fig. 1) to a single level and the same for the 3p sublevels. We have taken into account all the 
processes (collisional and radiative) which can populate or depopulate each of the levels. 

The resolution of the system of three equations gives us the ratios nL , with L=3d,4s,4p, being n Ca + the total 

n Ca + 

density of Ca + ions. The density and temperature we considered were 300cm -3 and 10000K, respectively. From this 
we deduced the relative population with respect the ground state. The results turn out to be: 

™ ~ io- 7 

and 

n ip 

~ U 

The negligible population of the upper 4p level prevents any contribution by cascade to the population of the 
3d level, therefore collisional excitation is the only important mechanism populating the metastable 3d level. This is 
consistent with the fact that we do not observe the triplet of Ca + (8498,8542,8662) (4p to 4s) in any EELR although 
it is observed in the broad line region of AGN where densities are higher by many order of magnitudes. 

a) Ionization of Ca+* by soft continuum photons. 

The soft UV counterpart of the ionizing continuum provides a source of ionization for both Ca + (IP: 11.9eV) and 
Ca + * (IP: 10.2eV). To estimate the photoionization rates, we will make the following approximaptions: 

1) At the fairly large depth in the cloud where the specie Ca + becomes abundant, we only need to consider photons 
with energies < 13.6eU, the ionization potential (IP) of H°, because photons just above this energy have already been 
absorbed and also because of the rapid decrease of the photoionization cross section with increasing energies. 

2) The ionizing continuum which reaches the PIZ is considered to be the soft UV counterpart of an ionizing PL of 
index —1.4 (but unattenuated since the opacity due to dust or trace elements below 13.6eV is relatively small). The 
continuum impinging the cloud is described by F v = F s v~ 1A . In number of photons this is F v jhv = F s //i^~ 2 - 4 . The 
constants F s /h will cancel out when taking the ratio n 3( i/L T 4 S . 

The atomic data was taken from Osterbrock (1987) for H° and from Shine & Linsky (1975) for Ca + and Ca+*. 
Tables 1, 2, 3 where we define a' v — a v * 10 18 and v' = z//10 16 show the repevant atomic data. The threshold ionizing 
frequency of H° is labelled vq. 

We estimate the quotient n 3( j/n4 S as follows: 

J v- 2A a v {Ca+*)dv 



4s / v-^ ai/ {Ca+)dv 

Vis 

If we define a' v — a v * 10 18 and v' = z^/10 16 , we have 

/ v'- 2A a' v {Ca+*)dv' 

ITw _ nsd^id 

ILu n is v' 

J v'- 2A a' v {Ga+)dv> 

= — t ( eq - 2 > 

The values of a' u (Ca + *) at v' , v' is and a' v (Ca + ) at v' a and v' is have been obtained by interpolation from (see Tables 2 
and 3). Using the rule of the rectangle to approximate a given integral: 

j f(x)dx ~ (xi - x )f(x ) + (x 2 - xi)f(xi) + ... 

+ (x n - X„-l)f(x„-l) 



n 4s 
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Table 1. H photoionization cross section 
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0.015 
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0.015 


0.001074 


3.2 


0.01 


0.0006133 


3.4 


0.01 


0.0005302 


3.6 


0.01 


0.0004622 



(eq.2) reduces to 

n 3d ^ n 3d 10.06 

II4 S 7l4 s 0.156 

If we take into account that ^ ~ 10~ 7 we obtain 

»4s 



(eq.3) 



^~6x 10- 6 

n 4s 

which means that the ionization of Ca+ into Ca ++ originates overwelmingly from the ground level 4s. 
b) Ionization of Ca+* by Lya photons. 

Lya photons have energies slightly higher than the treshold energy of level 3d of Ca + * and constitute undoudtedly 
an important source of ionization for this level. To compute the Lya emissivity, we assume the low density regime 
whereby ~ | of recombinations of H + lead to the emission of a Lya photon. Since Lya is a resonnant line of large 
line scattering opacity, we consider that the density of Lya photons within the Ca + region result from the production 
of Lya either locally or from deeper regions. The justification for this is that resonnant line photons generated from 
layers nearer the slab'surface would be reflected outward as a result of the increasing fraction of H° with depth (see 
Binette et al 1993b). Therefore the Lya flux potentially available to ionize level 3d is a small fraction 77 of the total 
Lya flux emitted by the cloud. This fraction is of order 0.2 corrresponding to the fraction of ionizing photons of H° 
not yet absorbed at the typical depth where the Ca + specie is abundant. 

We first compute the H^ a /H is ratio taking into account that resonnant scattering will increase the density 
of locally emitted line photons by a factor £ ~ 10 7 , the mean number of scatterings before escape. Such a high 
number of line scattering, however, characterizes only the locally produced Lya photons. This accumulation (or 
slowing down) effect which we want to estimate is only effective within a zone of optical depth of order a few. Let's 
take T scat ~ 10 = a scat n(H°) 5X where SX is the geometrical depth and n(H°) the local density of H°. Adopting 
a scat — 6 x 10~ 14 cm 2 (cf Appendix B of Binette et al 1993b), the column density of recombining H + we should consider 
in generating Lya is within a thickness N(H+) ~ N(H°) = n(H a ) 5X = 10/6 x 10~ 14 = 1.7 x 10 14 cm - 2 . 

Taking these considerations into account, the problem reduces to estimating the quotient 
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Table 2. Ca + 4s photoionization cross section 
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Table 3. Ca + * 3d photoionization cross section 
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2„ 



r£a Lya (Ca+*)N(H+) J v'- 2A a' p (H*)d v ' 

3d 



nLya 
3d n 3d 



n 4s n is v ' 

J v'- 2A a' v {Ca+)d V ' 



= ^\rta Lva {Ca+*)N{H+) I -± (eqA) 
ri4, s 6 I2 

where a_L yQ (Ca + *) is the photoionizatiom cross section of Ca+* at the Lya energy (~ 6 x 10~ 18 cm 2 ). Taking again 
n^d/n^s <~ 10~ 7 , we obtain 

=f«L_ ~ 1.36 x 1(T 4 ^ (eg.5) 
His h 

Using again the rule of the rectangle we obtain 



n 



3d - 1.5 x 1(T 2 . 



n 4s 

The effect of slowing down of resonnant line photons is interesting but appears insufficient as it envolves too small 
a fraction of the Lya photons generated within the nebula. Let's now estimate the importance of all the Lya photons 
generated within the deeper zones which, after having scattered far enough in frequency to escape, must still cross the 
Ca + zone (without appreciable scattering in that zone) . The effect on the photoionization of level 3d is given in this 
case by 

00 

ha' Lya (Ca+*) f v>-™dv> 

3a- = ^— * 

n 4s n is " 

J v'- 2A a' v {Ca+)dv' 



n^d 2 , , +*^4 , _a 
n 4s 3 y h 

(Note that a' is used here for Ca + *). Using the rule of the rectangle we obtain 

2.43 x 10" 6 

In summary, Lya emission and/or trapping are insufficient to ionize Ca + . As in the case of soft continuum photons, 
this is basically the result of the extremely small population characterizing the excited level 3d. 
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